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Some possibilities to test the Higgs sector of the Minimal Supersymmetric Standard Model 
(MSSM) at the LHC and future e + e _ colliders are discussed. This includes precision coupling 
strength measurements, the search for additional Higgs bosons as well as their decay to supersym¬ 
metric particles. 


I. INTRODUCTION 

The ATLAS and CMS experiments at CERN have discovered a new boson with a mass around 125 GeV EL@- 
Within the present experimental uncertainties this new boson behaves like the Higgs boson of the Standard 
Model (SM) [3|. However, the newly discovered particle can also be interpreted as the Higgs boson of extended 
models, where the Minimal Supersymmetric Standard Model (MSSM) [4] is a prime candidate. The Higgs sector 
of the MSSM with two scalar doublets accommodates five physical Higgs bosons. In lowest order these are the 
light and heavy CP-even h and H , the CP-odd A, and the charged Higgs bosons H ± . It can be expressed (at 
lowest order) in terms of the gauge couplings, the mass of the CP-odd Higgs boson, Ma, and tan/3 = V 2 /V 1 , 
the ratio of the two vacuum expectation values. All other masses and mixing angles can therefore be predicted. 
Higher-order contributions can give large corrections to the tree-level relations @]. An upper bound for the 
mass of the lightest MSSM Higgs boson of Mh 135 GeV had been obtained [?], in perfect agreement with the 
observed value of 


M e * p = 125.09 ± 0.24 GeV , (1) 

as evaluated by the combination of ATLAS and CMS measurements f8)]. 

We will review a few ways to test the MSSM Higgs sector at the LHC and beyond. We first briefly discuss 
the precision prediction for the lightest CP-even Higgs boson mass in the MSSM. We review the status and the 
prospects of the Higgs coupling strength analyses at the LHC and the ILC. Finally we discuss where additional 
MSSM Higgs bosons could be discovered and review some precision calculations for their (potential) decay to 
supersymmetric (SUSY) particles. 


II. THE LIGHTEST HIGGS BOSON MASS AS A PRECISION OBSERVABLE 


In the MSSM the mass of the light CP-even Higgs boson, Mh, can directly be predicted from the other 
parameters of the model. The accuracy of this prediction should at least match the one of the experimental 
result. The measured Higgs-boson mass value, Eq. ©, has already reached the level of a precision observable 
with an experimental accuracy of about 250 MeV. Consequently, it plays an important role in the context of 
testing the MSSM Higgs sector. 


The status of higher-order corrections to Mh is quite advanced, see Refs. M 10] for the calculations of the 
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full one-loop level. At the two-loop level |ll|, |L2j, |l 
particular the 0(a t a s ) and 0(a 2 ) contributions (a t = h 2 /(Ai r), h t being the top-quark Yukawa coupling) to the 
self-energies - evaluated in the Feynman-diagrammatic (FD) as well as in the effective potential (EP) method 
- as well as the 0(aba s ), 0(atab) and 0(a 2 ) contributions - evaluated in the EP approach - are known for 
vanishing external momenta. An evaluation of the momentum dependence at the two-loop level in a pure DR 
calculation was presented in Ref. [27J . A (nearly) full two-loop EP calculation, including even the leading three- 
loop corrections, has also been published H^]. However, the calculation presented in Ref. [28:] is not publicly 
available as a computer code for Higgs-mass calculations. Subsequently, another lea ding three-loop calculation 
of 0(ata 2 s ), depending on the various SUSY mass hierarchies, has been performed [29j ] . resulting in the code 
H3m which adds the three-loop corrections to the FeynHiggs [H, HQ;, HH, [32] result. Recently, a combination 
of the full one-loop result, supplemented with leading and subleading two-loop corrections evaluated in the 
Feynman-diagrammatic/effective potential method and a resummation of the leading and subleading logarithmic 
corrections from the scalar-top sector has been published [32j in the latest version of the code FeynHiggs 

QH \M-\M- 

More recently the calculation of the momentum dependent two-loop QCD corrections to Mh have been 
presented [33]. (From a technical point of view we have calculated the momentum dependent two-loop self- 
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energy diagrams numerically using the program SecDec [H,[HiEll-) Su bseq uently, in Ref. [33] this calculation 
was repeated (differences of the two calculations are discussed in Ref. :38j), where also a calculation of the 
two-loop corrections of 0(aa s ) were presented. The results of Ref. [33 are publicly available in the code 
FeynHiggs. 

The remaining theoretical uncertainty in the calculation of Mh, from unknown higher-order corrections, had 
been estimated to be up to 3 GeV, depending on the parameter region. Recent improvements have potentially 
lead to a somewhat smaller estimate of up to ~ 2 GeV [32j, [3f| for not too large SUSY mass scales. However, 
a careful re-analysis of this uncertainty for lower and heavier SUSY mass scales is in order. As the accuracy of 
the Mh prediction should at least match the one of the experimental result, further sub-dominant higher-order 
corrections have to be included in the Higgs-boson mass predictions [I(| • 


III. COUPLING STRENGTH ANALYSIS AT THE LHC AND BEYOND 

Testing the coupling strengths of the discovered Higgs boson could yield hints towards an extended Higgs 
sector, where the MSSM Higgs sector makes clear predictions for possible deviations. In order to test the 
compatibility of the predictions for the SM Higgs boson with the (2012) experimental data, the LHC Higgs Cross 
Section Working Group proposed several benchmark scenarios for “coupling scale factors” [4l|, [42| (see Ref. [43j] 
for a recent review on Higgs coupling extractions). Effectively, the predicted SM Higgs cross sections and partial 
decay widths are dressed with scale factors (and = 1 corresponds to the SM). Several assumptions are 
made for this re-framework: there is only one state at 125 GeV responsible for the signal, the coupling structure 
is the same as for the SM Higgs (i.e. it is a CP-even scalar), and the zero width approximation is assumed to be 
valid, allowing for a clear separation and simple handling of production and decay of the Higgs particle. The 
most relevant coupling strength modifiers are re t , re&, re T , kw , «z, re 7 , re s , .... 

One limitation at the LHC (but not at the ILC) is the fact that the total width cannot be determined 
experimentally without additional theory assumptions. In the absence of a total width measurement only ratios 
of re’s can be determined from experimental data. An assumption often made is nw,z < 1 33]- A recent analysis 
from CMS using the Higgs decays to ZZ far off-shell yielded an upper limit on the total width about four times 
larger than the SM width (45[. However, here the assumption of the equality of on-shell and off-shell couplings 
of the Higgs boson plays a crucial role. It was pointed out that this equality is violated in particular in the 
presence of new physics in the Higgs sector [43, S3] ■ 

In the left plot of Fig. [T] we compare the results estimated for the HL-LHC (with Sab” 1 and an assumed 
improvement of 50% in the theoretical uncertainties) with the various stages of the ILC under the theory 
assumption kw,z < 1 Ssj- This most general fit includes kw,z for the gauge bosons, n U} d,i for up-type quarks, 
down-type quarks and charged leptons, respectively, as well as re 7 and re g for the loop-induced couplings of the 
Higgs to photons and gluons. Also the (possibly invisible) branching ratio of the Higgs boson to new physics 
(BR (H —> NP)) is included. One can observe that the HL-LHC and the ILC 250 yield comparable results. 
However, going to higher ILC energies, yields substantially higher precisions in the fit for the coupling scale 
factors. In the final stage of the ILC (ILC 1000 LumiUp), precisions at the per-mille level in nw,z are possible. 
The 1 — 2% range is reached for all other re’s. The branching ratio to new physics can be restricted to the 
per-mille level. 

Using ILC data the theory assumption Kw,z < 1 can be dropped, since the “Z-recoil method” (see Ref. [49| 
and references therein) allows for a model independent determination of the HZZ coupling. The corresponding 
results are shown in the right plot of Fig. |TJ where the HL-LHC results are combined with the various stages 
of the ILC. The results from the HL-LHC alone continue to very large values of the re’s, since the fit cannot 
be done without theory assumptions. Including the ILC measurements (where the first line corresponds to the 
inclusion of only the cr^f 1 measurement at the ILC) yields a converging fit. In the final ILC stage nw,z are 
determined to better than one per-cent, whereas the other coupling scale factors are obtained in the 1 — 2% 
range. The branching ratio to new physics is restricted to be smaller than one per-cent. This opens up the 
possibility to observe MSSM induced deviations in the Higgs boson couplings, provided that the overall Higgs 
mass scale, Ma , is not too large. 


IV. THE SEARCH FOR ADDITIONAL HIGGS BOSONS 

Many investigations have been performed analyzing the agreement of the MSSM with a Higgs boson at 
~ 125 GeV. In a first step only the mass information can be used to test the model, while in a second step also 
the rate information of the various Higgs search channels can be taken into account (see the previous section). 
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FIG. 1: Fit to the coupling scale factors with (left) and without the theory assumption of Kw,z < 1 


Here we briefly review some results in two of the new benchmark scenarios (5Q), devised for the search for 
heavy MSSM Higgs bosons. In the left plot of Fig. [2] the m™ ax scenario is shown. The red area is excluded 
by LHC searches for the heavy MSSM Higgs bosons, the blue area is excluded by LEP Higgs searches, and 
the light shaded red area is excluded by LHC searches for a SM-likc Higgs boson. The bounds have been 
obtained with HiggsBounds [51] (where an extensive list of original references can be found). The green area 
yields Mh = 125 ± 3 GeV, i.e. the region allowed by the experimental data, taking into account the theoretical 
uncertainty in the Mh calculation as discussed above. Since the m™ ax scenario maximizes the light CP-even 
Higgs boson mass it is possible to extract lower (one parameter) limits on Ma and tan /? from the edges of the 
green band. By choosing the parameters entering via radiative corrections such that those corrections yield a 
maximum upward shift to Mh, the lower bounds on Ma and tan/3 that can be obtained are general in the sense 
that they (approximately) hold for any values of the other parameters. To address the (small) residual Msusy 
dependence (-Msusy denotes the average scalar top mass scale) of the lower bounds on Ma and tan/?, limits 
have been extracted for the three different values Msusy = 10.5,1,2} TeV, see Tab. [I] [52|. For comparison 
also the previous limits derived from the LEP Higgs searches [531 are shown, i.e. before the incorporation of the 
Higgs discovery at the LHC. The bounds on Ma translate directly into lower limits on M H ±, which are also 
given in the table. More recent experimental Higgs exclusion bounds shift these limits to even higher values, 
see the left plot in Fig. [2] Consequently, the experimental result of Mh ~ 125 ± 3 GeV requires M H ± > m t 
with important consequences for the charged Higgs boson phenomenology. 

In the right plot of Fig. [2] we show the rn™ od+ scenario that differs from the m“ ax scenario in the choice of 
X t (the off-diagonal entry in the scalar top mass matrix). While in the m™ ax scenario X t /Msusy = +2 had 
been chosen to maximize Mh, in the m™ od+ scenario X t /Msusy = +1-5 is used to yield a “good” Mh value 
over the nearly the entire Ma ~tan /? plane, which is visible as the extended green region. 


Msusy (GeV) 

Limits without Mh ~ 125 GeV 
tan/3 M a (GeV) M H± (GeV) 

Limits with Mh ~ 125 GeV 
tan/3 M a (GeV) M H± (GeV) 

500 

2.7 

95 

123 

4.5 

140 

161 

1000 

2.2 

95 

123 

3.2 

133 

155 

2000 

2.0 

95 

123 

2.9 

130 

152 


TABLE I: Lower limits on the MSSM Higgs sector tree-level parameters Ma ( M H ±) and tan /? obtained with and without 
the assumed Higgs signal of Mh ~ 125.5 GeV. The mass limits have been rounded to 1 GeV [s3 ]. 
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FIG. 2: Myt-tan/3 plane in the m™ ax scenario (left) and in the m™ od+ scenario (right) [5C|. The green shaded area yields 
Mh ~ 125 ± 3 GeV, the red area at high tan /3 is excluded by LHC heavy MSSM Higgs boson searches, the blue area is 
excluded by LEP Higgs searches, and the red strip at low tan/3 is excluded by the LHC SM Higgs searches. 


V. PRECISION PREDICTIONS FOR THE DECAY OF HIGGS BOSONS TO SUSY PARTICLES 

Depending on the scale of Higgs and SUSY masses the main decay channels of the additional Higgs bosons 
could go to SUSY particles, which is demonstrated in Fig. [3] [5(j. The branching ratios for the decay of H 
and A into charginos and neutralinos may become large at small or moderate values of tan (3. In Fig. [3] we 
show the nr™ od+ (left) and m™ od ~ (right) scenarios [501, where the masses of the charginos and neutralinos 
are 0(200 GeV). The excluded regions from the Higgs searches at LEP and the LHC are as before. The color 
coding for the allowed region of the parameter space indicates the average value of the branching ratios for 
the decay of H and A into charginos and neutralinos (summed over all contributing final states). One can see 
from the plots that as a consequence of the relatively low values of the chargino/neutralino masses in these 
benchmark scenarios the decays of H and A into charginos and neutralinos are kinematically open essentially 
in the whole allowed parameter space of the scenario, with the exception of a small region with rather small 
Ma- The branching ratios for the decays of H and A into charginos and neutralinos reach values in excess of 
70% for small and moderate values of tan /3. Including these channels into the searches for heavy MSSM Higgs 
bosons could potentially allow to discover new Higgs bosons and SUSY particles at the same time. 

Recently, full one-loop calculations for the decays of Higgs bosons to scalar fermions [54| and into 
charginos/neutralinos [55| in the MSSM with complex parameters (cMSSM) have become available. In Fig. 0] 
we show the results for the decay hi —> X 2 X 2 (* = 2,3), where details and parameter settings can be found in 
Ref. [5l]. /12 and /13 are the two neutral heavy Higgs bosons in the cMSSM, corresponding to H and A in the 
real case. In the left plot of Fig. 2] the decay widths at the tree- and at the full one-loop level are shown as 
a function of M H ± (at M H ± ~ 1000 GeV and M H ± ~ 1520 GeV a mass crossing of hi and /13 takes place, 
see Ref. [55| for details). It can be seen that the (in this case purely electroweak) one-loop correction can 
change the decay width by up to 20%. In the right plot of Fig. [4] the decay widths are shown as a function of 
¥>Md the phase of the 0(1) gaugino soft SUSY-breaking parameter. Changing the phase can lead to effects of 
up to 50%, while the one-loop corrections again can be as large as 20%. These examples show that complex 
parameters and the full one-loop corrections should be taken into account for the interpretation of the searches 
for charginos/neutralinos as well as for any future precision analyses of those decays. 
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FIG. 3: JWA-tan/? plane in the m r f od+ scenario (left) and the m™ od_ scenario (right) [50[. he exclusion regions are shown 
as in Fig. [2] while the color coding in the allowed region indicates the average total branching ratio of H and A into 
charginos and neutralinos. 
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FIG. 4: T(hi —>• X 2 X 2 )• Tree-level and full one-loop corrected partial widths for hi —1 X 2 X 2 (* = 2,3) are shown [55il . 
The left plot shows the partial decay width with M H ± varied. The right plot shows the complex phase ipM t varied (see 
text). 
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